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Abstract Recently, we have derived a two-nucleon potential and consistent nuclear electromagnetic currents 
in chiral effective field theory with pious and nucleons as explicit degrees of freedom. The calculation of 
the currents has been carried out to include N^LO corrections, consisting of two-pion exchange and contact 
contributions. The latter involve unknown low-energy constants (LECs), some of which have been fixed by 
fitting the np S- and P-wave phase shifts up to 100 MeV lab energies. The remaining LECs entering the current 
operator are determined so as to reproduce the experimental deuteron and trinucleon magnetic moments, as 
well as the np cross section. This electromagnetic current operator is utilized to study the nd and n"^He 
radiative captures at thermal neutron energies. Here we discuss our results stressing on the important role 
played by the LECs in reproducing the experimental data. 
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Quantum chromodynamics (QCD) is the under- 
lying theory of the strong interaction. On this basis, 
interactions among the relevant degrees of freedom of 
nuclear physics, such as pions, nucleons, and delta- 
isobars, are completely determined by the quark and 
gluon dynamic. At low energies though, the strong 
coupling constant becomes too large to allow for ap- 
plication of perturbative techniques to solve QCD. 
Consequently, we are still far from a quantitative un- 
derstanding of the low-energy physics by ab initio 
calculations from QCD. Chiral effective field theory 
(xEFT) exploits the symmetries exhibited by QCD 
in the low-energy regime, in particular chiral symme- 
try, to constrain the form of the interactions of the 
pions among themselves and with the other degrees 
of freedom'^'. The pion couples by powers of its mo- 
mentum Q and the Lagrangians describing these in- 
teractions can be expanded in powers of Q/A^, where 
A^ ~ 1 GeV represents the chiral-symmetry breaking 
scale and characterizes the convergence of the expan- 
sion. The effectiveness of the theory is then confined 
to kinematic regions where the constraint Q ^ A^ is 



realized. The unknown coefficients of the chiral ex- 
pansion, i.e. the low energy constant (LECs), need to 
be fixed by comparison with the experimental data. 
^EFT provides an expansion of the Lagrangians in 
powers of a small momentum as opposed to an ex- 
pansion in the strong coupling constant, restoring de 
facto the applicability of perturbative techniques also 
in the low-energy regime. Due to the chiral expansion 
it is possible, in principle, to evaluate an observable 
to any degree of desired accuracy and to know a pri- 
ori the hierarchy of interactions contributing to the 
low energy process under study. 

Since the pioneering work of Weinberg'^', this cal- 
culational scheme has been widely utilized in nuclear 
physics and nuclear xEFT has developed into an in- 
tense field of research. Nuclear two- and three-body 
interactions^^' , as well as interactions of electroweak 
probes with nuclei'^' ^ have been studied within the 
xEFT approach. 

Recently, we have derived the nuclear electromag- 
netic (EM) currents in xEFT'^'-^, retaining, as de- 
grees of freedom, pions and nucleons. The calcula- 
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tion has been carried out in time-ordered perturba- 
tion theory'^ with non-rclativistic Hamiltonians de- 
rived from the chiral Lagrangians of Refs. 0, S, 
The strong and electromagnetic interaction Hamilto- 
nians required to evaluate the EM current operator 
up to N^LO accuracy — that is eQ in the chiral ex- 
pansion, Q denoting the low momentum scale, and e 
being the electric charge — are listed in Ref. 0, 01 ■ 

In Fig. [1] we show the contributions to the current 
operator up to N^LO (eQ°). The LO {eQ-^) term 
is given by a one-body contribution, consisting of the 
standard convection and spin-magnetization nucleon 
currents, while pion-exchange currents occur at NLO 
(eQ-i). The N^LO term is due to {Q/Mf relativistic 
corrections — where M denotes the nucleon mass — to 
the LO one-body current. 

In Fig. [5] we list the N'^LO contributions, which 
can be separated into three classes: i) one-loop 
two-pion exchange terms, represented by diagrams 
(a)-(i); ii) tree- level term involving the nuclear- 
electromagnetic Hamiltonian of order e at the ver- 
tex illustrated by a full circle in diagram (j); and iii) 
contact currents of minimal and non-minimal nature, 
illustrated by diagram (k). 

The last two contributions involve unknown 
LECs. In particular, the tree-level current of the type 
shown in panel (j), depends on three LECs, two of 
them multiply isovector structures and the remain- 
ing one multiplies an isoscalar structure. Incidentally, 
the isovector part of this tree-level current has the 



same structure as the current involving the excita- 
tion of a delta-isobar'*''. This resonance saturation 
argument is exploited to infer the ratio between the 
two LECs multiplying the isovector terms in the cur- 
rent of diagram (j) (see below). Contact currents of 
non- minimal character, panel (k) in Fig.[Tl depend on 
two additional unknown LECs, multiplying respec- 
tively an isoscalar and an isovector structure, while 
those obtained via minimal substitution are expressed 
in terms of LECs entering the contact two-nucleon 
chiral potential of order (or N^LO) ^ . The two- 
nucleon potential has been derived in Ref. Q up to 
N^LO and these LECs have been fixed by fitting the 
np S- and P-wave phase shifts up to 100 MeV lab- 
oratory energies'^. Thus total number of unknown 
LECs to be determined is reduced to four. 



LO : 

NLO : eQ-1 
N^LO : eQ° 




Fig. 1. Diagrams illustrating one- and two- 
body currents up to N^LO (eQ°). Nucle- 
ons, pious, and photons are denoted by solid, 
dashed, and wavy lines, respectively. The 
square represents the relativistic correction to 
the LO one-body current. Only one among 
the possible time orderings is shown for the 
NLO diagrams. 
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Fig. 2. Diagrams illustrating two-body currents entering at N"^LO (eQ), notation as in Fig. [T] Only one 
among the possible time orderings is shown for diagrams (a)-(j). 



An important aspect of the derivation of the EM 
currents (and two-nucleon potential) is to retain both 
irreducible diagrams and recoil-corrected reducible 
onesl^. The latter arise from expanding the energy 
denominators (in reducible diagrams) in powers of 
nucleon kinetic energy differences to pion energies 
(these ratios are of oder Q). Partial cancellations 
occur between the irreducible and recoil-corrected re- 
ducible contributions both at N^LO and N^'LO'®. We 
also note that this approach leads to N^LO EM cur- 
rents that satisfy the continuity equation with the 
corresponding N^LO two-body potential"''. The ex- 
pressions for the two-pion-exchange N^LO currents 



in panels (a)-(i) of Fig.[2]are in agreement with those 
obtained by KoUing et al. in Ref. [l3| by the method 
of the unitary transformations. However, they are 
different from those derived by Park et al. in Ref. Q 
in covariant perturbation theory, since these authors 
include irreducible contributions only. 

We now present a study of the nd and n^He ra- 
diative capture at thermal neutron energies within 
the hybrid approach, where the EM xEFT current 
operator described above is used to evaluate transi- 
tion matrix elements between nuclear wave functions 
obtained with realistic Hamiltonian with two- and 
three-body potentials. In order to study the model 
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dependence of the calculated observables, we use two 
different combinations of two- and three-body po- 
tential, namely the Argonne Vig^ with the Urbana- 
IX^ three-nucleon potential (AV18/UIX), and the 
N^LOtiS and N^LO'"''' chiral two- and three-nuclcon 
potentials (N3LO/N2LO). We study the sensitivity 
of the observables to variations of the cutoff A, in- 
troduced to regularize the EM current operator via 
the momentum cutoff CA(fc) = exp(— fc'*/A'*). In our 
study, A varies from 500 to 700 MeV which corre- 
sponds to "removing" short-range physics at distance 
scales less l/(3m^). 
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Fig. 3. Cumulative LO, NLO, N^LO, and 
N'^LO(S-L) contributions for the deuteron and 
trinucleon isoscalar and isovector magnetic 
moments, and np radiative capture. 



Out of the four unknown LECs entering the EM 
current operator, two multiply isoscalar structures 
and two multiply isovector operator structures. We 
fix these LECs by reproducing the experimental val- 
ues of two isoscalar observables, i.e. the deuteron 
[/id] and the isoscalar [/i^(^He/^H)] combination of 
the trinucleon magnetic moments, and two isovector 
observables, i. e. the isovector [/^^(^He/^H)] combi- 
nation of the trinucleon magnetic moments and the 
np cross section [cr^p] at thermal neutron energies. 
The results are shown in Fig. [3] where the cumulative 
contributions at LO, NLO, N^LO, and N3L0(S-L) are 
represented. The cumulative contribution N^LO(S-L) 
is given by the terms up to N^LO plus the N^LO con- 
tributions associated with pion loops (represented in 
panels (a)-(i) of Fig. (2), which depend on the (known) 
nucleon axial coupling constant, pion decay ampli- 
tude, and pion mass, as well as with contact currents, 
which depend on the LECs obtained from the fits to 
the np phase shifts. 

The LECs entering the complete current, denoted 
in what follows as N''LO(LECs), are fixed, for each 
value of the cutoff A, so as to reproduce the experi- 
mental values which in Fig. [3] are represented by the 



black band, including experimental errors. The sen- 
sitivity of the results to the two Hamiltonian models 
utihzed (AV18/UIX and the N3LO/N2LO) is repre- 
sented by the thickness of the color bands. We note 
that the sign of the N^LO and N^LO(S-L) contribu- 
tions is opposite to that of the LO and NLO contribu- 
tions. This increases the discrepancy between theory 
and experiment. 




500 600 700 
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Fig. 4. Cumulative LO, NLO, N^LO, N^LO(S- 
L), and N^LO(LECs) contributions to the nd 
(°"nd) ^^'^ n'^He (o"^3jjj,) cross sections (right 
and left top panel respectively), and circular 
polarization factor Rc. 



Having fixed all the LECs, we are left with a com- 
pletely determined EM current operator which can 
now be used to make predictions for the n(d,7)^H 
and 7T,(^He,7)*He reactions' cross sections — denoted 
as (7^^ and cr^sjj^, respectively — and the circular po- 
larization factor associated with the capture of po- 
larized neutrons on deuterons. In this calculation we 
have used the AV18/UIX (N3LO/N2LO) combina- 
tion of two- and three-nucleon potentials for the A —3 
(yl=4) processes; calculations with the N3LO/N2LO 
(AVI8/UIX) potential models are in progress. The 
predictions are represented in Fig. [3] along with the 
experimental data, shown in black, which are from 
Ref. ^15^ for nd and Ref. 16 1 for n^He. The com- 
plete N^LO(LECs) current is shown in Fig. H] by the 
orange lines. The calculated nd cross section is in 
excellent agreement with the measured value and is 
weakly dependent on the cutoff. The cross section 
for the n^He reaction undergoes a 5% variation when 
the cutoff changes from 500 to 700 MeV, but is still 



No. X 



L. Girlanda et al: Electromagnetic processes in a xEFT framework 



4 



within the experimental error band. These reactions 
are known to be dominated by many-body compo- 
nents of the current operator, which provide most of 
the calculated cross section^^S. This trend is con- 
firmed here: the LO contribution to the cross sec- 
tions is highly suppressed, and provides only about 
46% (18%) of the total calculated nd (n^He) value. 
What is more interesting though, is the large contri- 
bution associated with the N^LO(LECs) currents in 
both these reactions. These currents are crucial for 
bringing theory into agreement with experiment. 

We are presently in the process of extending these 
hybrid studies to different realistic Hamiltonian mod- 
els, with the goal of quantifying the sensitivity of the 
cross sections to the wave functions employed in the 
calculations. Obviously, our ultimate objective is to 
perform a fully consistent xEFT calculation, using 
the N^LO potential derived in Ref. [7], along with 



the EM currents we presented here. In Ref. [7| we 
show the deuteron wave functions obtained with the 
N^LO chiral potential and compare them with those 
corresponding to the AV18. The two sets of wave 
functions display a different behavior at short range, 
in particular the N^LO D-wave component is signifi- 
cantly smaller than the AVIS. From this perspective, 
it will be interesting to establish whether these chiral 
potential and currents lead to a satisfactory descrip- 
tion of the nd and n^He captures. 
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by the U.S. Department of Energy, Office of Nu- 
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